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ABSTRACT 

The enhanced star formation in the inner 100 pc of the Galaxy launches a superwind, reaching ~ 1600 km s" 1 
for M82-like parameters. The ram pressure of the wind is very low compared to more powerful starburst winds. 
I show that halo gas stops the wind a few kpc from the Galactic Center. I argue that the termination shock 
accelerates cosmic rays, and that the resulting Inverse Compton 7-rays are visible as the Fermi Bubbles. The 
Bubbles are thus wind bubbles, which the starburst can inflate within 10 Myr. They remain in steady state as 
long as the starburst lasts. The shock may accelerate PeV electrons and ultra-high energy protons. The Bubbles 
may be analogs of galactic wind termination shocks in the intergalactic medium. I discuss the advantages and 
problems of this model. 

Subject headings: Galaxy: center — galaxies: starburst — ISM: jets and outflows — gamma rays: galaxies — 
cosmic rays 
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1. INTRODUCTION 

Powerful outflows are ever present in starburst regions 
(Ssfr <^ 0.1 M0 yr" 1 kpc" ), blasting out at speeds 3> 
100 k m s" 1 dChevalier & Clegg||1985l( CC85):lHeckman et all 
Il990t IStrickland & Heckmari 120091 Winds erupt in 
these intense regions when superbubbles and supernova 
remnants overlap as the star-formation density increases 
dMcKee & O striker 1977). The hot plasma fills starbursts and 
blows out in a sound-crossing time, continuously replenished 
by massive stars. Even in normal galaxies like the Milky 
Way, cosmic rays (CRs) can drive winds once they diffuse fa r 
enough from the galactic plane dBreitschwerdt et aT1ll991l) . 
But what happens to the wind after it escapes these central 
regions? 

In the CC85 model, the wind rapidly accelerates as it 
passes the sonic point, reaching an asymptotic speed Voo ~ 
1600 km s if supernova energy thermaliz es efficiently and 
the m ass loading is small as in M82 ( Strickland & Heckman 
120091) . As the wind expands adiabatically out beyond the con- 
fines of the starburst, particle random motions in the wind 
slow down. This appl ies both to the w ind plasma itself and 
CRs within the winds (IVolket al.ll 19961) . The random kinetic 
energy of the wind is converted into bulk kinetic energy, push- 
ing the wind plasma like a piston. 

The wind usually escapes the galaxy, growing ever more di- 
lute. Eventually, it becomes so rarefied that the ram pressure 
equals whatever external pressure there is. At this location, a 
termination shock stops the wind. These shocks are theorized 
to host energetic phenomena, such as cosmic ray acceleration. 
However, the termination shock of starbursts lie in the inter- 
galactic medium. 

A prototype for starburst regions and outflows exists in our 
Galaxy, in the form of the Galactic Center Central Molecular 
Zone (CMZ). Lying within 100 pc o f Sgr A*, the region is 
filled with gas (Molinari et al.ll201lh . a large star- formation 
rate d ensity (Esfr ~ 2 M Q yr" 1 kpc" 2 : lYusef-Z adeh et al. 
120091) . str ong magnetic fields fcrockeret al.ll2010l) , and TeV 
emission dAharonian et al. 2006]). Evidence for an outflow in- 
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elud es the abnormally wea k nonthermal radio and 7-ray emis- 
sion dCrocker et al.ll201 la[). suggesting that CRs are advected 
away from the CMZ dCrocker et al.ll20lTbb . 

The Fermi Bubbles are mo re evidence for powerful phe- 
nomena in the Galact i c Center (Finkbeiner 2004; Dobler et al. 
l2010t lSuetalJl2O!0t lAde et all 120 121) . These bilobal struc- 
tures, visible in 7-rays and radio, extend nearly 10 kpc from 
the Galactic Plane and emanate from the Center. Soft X-ray 
emission is suggestive of a shock at the Bubble edges dSofud 
l2000HBland-Hawthorn & Cohenll200H [Su et alj|2010b . Most 
theoretical work interprets these bubble s as the result of an 
outburst from Sgr A* (Cheng et al. 201 1; Zubovas et a 0201 It 
IGuo & Mathewsll2012l) . In these models. CR electrons (e") 
shine in radio by synchrotron emission and in 7-rays by In- 
verse Compton (IC) emission. But, since these losses are 
rapid, either the <?" are transported rapidly or they must be 
a ccelerated in place. 

iCrocker & Aharonianl(l201 11) proposed an alternative expla- 
nation, in which the CMZ powers the Bubbles. In this model, 
the CMZ has been accelerating CRs for the entire Galaxy's 
history. The wind advects CR protons a few kpc, where they 
accumulate for ten billion years. The Bubbles' 7-ray emis- 
sion is pionic; secondary e ± made by pionic interactions emit 
the radio waves. Since the gas density is low in the Bubbles 
(< 0.01 cm" 3 ), the CRs must accumulate for ten billion years 
in order to make enough 7-rays. This model evades the se- 
vere adiabatic losses in a CC85-style wind by postulating a 
slow wind that becomes compress ed in the Bubbles . Lepto nic 
starburst models are possible too: Biermann et al. (12010b ar- 
gued that high energy e ± rapidly diffuse from the Center's 
starburst and are advected into the Bubbles, where they emit 
IC. 

In this Letter, I propose that the Fermi Bubbles are the ter- 
mination shocks of the Center's starburst wind (Figure [T). 
Because the CMZ wind is less dense than prototypical star- 
bursts' winds, its ram pressure is quite low. Thus, the wind is 
stopped by the Galactic halo. These shocks accelerate CRs. 
I suggest that primary <?" are the sources of the nonthermal 
emission. Thus, the Fermi Bubbles are a nearby analog of the 
inter galactic termination shocks of more pow erful starbursts 
(c.f.. Uokipii & Morfill| [T987; Vol k & Zirakas hvili 2004]). 




FIG. 1 . — Sketch of the CMZ and its termination shock with respect to the 
Galaxy. CR e~ and p + are accelerated at the shock; the emission from the e~ 
appears to us as the radio and 7-ray bubbles. 

2. PROPERTIES OF THE CMZ WIND AND ITS TERMINATION 
SHOCK 

2.1. Basic wind properties 

My discussion is informe d by the CC85 theory of star- 
burst winds, as amended by IStrickland & Heckmanl (12009b 
to apply to disk geometries. A hot plasma fills the wind 
volume, with properties determined by energy and mass in- 
jection. The CMZ is a disk with radius Rcmz = 112 pc 
and height /jc mz = 42 pc, forming stars at a rate SFR « 
0.1 M Q yr" 1 (ICrocker et all 1201 lbT) . Mass is injected into 
the wind at a rate M = 0.117/3 x SFR from stellar winds 
and explosions, and these phenomena heat the wind at a rate 
E = 2.54 x 10 41 e,hemJSFR/(M p yr^lerg sec" 1 fl I take = 2 
and etherm = 0.75, as for M82 (Strickland & Heckman 2009). 
Then the central density of the wind is given by 



0.592 JM 3 /E 
n r = — y — = 0.037 cm -3 



d2 /1 I 2/icmz 
K CMZ^ i+ Rcmz 



)pm H 



SFR 



O.lM yr" 



and the central temperature is 



r c= (7-iW? =3-7xl0 7 K 
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for a mean molecular weight p = 0.59 and adiabatic index 7 = 
5/3. 
In the CC85 model, the wind flows out of the starburst at 



a speed Voo = \J2E/M = 1600 km s '. The escape speed 

from the CMZ is only 800 - 1000 km s" 1 ( iBreitschwerdt et al.l 
U99lHLaunhardt et al.ll2002l) . so the wind escapesQ Although 
CC85 assume spherical symmetry, actual starburst winds are 



3 I have converted the values in Strickland & Heckman ( 2009) to a Salpeter 
IMF running from 0. 1 Mq to 100 Mq . 

4 Carretti et al. (2013) argue that the radio "spurs" trace CMZ features if 
the outflow has a speed of 1000 km s . 



observed to flow in cones. I therefore assume the outflow 
cone fills a solid angle CI, The density of the outflow is 

p = M/iCtVoor 2 ), and the ram pressure is P mm = pv^/2. 

2.2. The Termination Shock 

The termination shock occurs when the ram pressure is 
equal to the external pressure. Supposing the external pres- 
sure is entirely thermal, P e xt = n ext kBT ext , the termination shock 
distance is 



R,= 



Mv n 



2fln ext k B T ey 



(3) 



The state of the gas several kpc into the Galactic halo, and thus 
the location of the CMZ termination shock, is poorly known. I 
take P ex t/kB = 1000 K cm" 3 , consistent with loo se constraints 
from Galactic absorption line features (e.g., ISavage et al.l 
12001 lYao & Wandl200lt IHsu et al.l|201 11) and hydrodynam - 
ics (e.g.. iSpitzeril 19561: iWolfire et al.lll995l iFang et al.ll2013l) . 
R, depends weakly on the external pressure. 
For the CMZ wind, 



R, = 5A kpc 
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roughly the size of the Fermi Bubbles. For SFRs of 3 to 200 
M yr" 1 (e.g., NGC 253, M82, or Arp 220), the termination 
shock is instead firmly in the realm of circumgalactic or even 
intergalactic distances. 

The termination shock is a steady-state feature, present 
even if the starburst wind has always been active. The wind 
just inside the shock is in pressure equilibrium with the gas 
outside. In contrast, many models of the Fermi Bubbles pro- 
pose that the shock is transient, resulting from a recent burst 
of energy injection, and the shock is driven outwards by an 
overpressure. 

By the time the wind reaches these distances, it is very 
cold. Outside of the starburst, where there is no energy in- 
jection into the wind, the temperature of the wind scales as 
T oc p 7 " 1 oc p 2 / 3 ; likewise the energy density of relativistic 
CRs scales as p 1 / 3 . After scaling to the central density and 
temperature, 



T(R)^ 1.7 x 10 5 K 



R 



Ucr(R) 
Ucr(0) 
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The adiabatic cooling time {2,/A)R/v OQ w 2.8 Myr (R/b kpc) 
is comparable to the radiative co oling time for high (^ 3Z ) 
metallicity gas at these distances (ISufherland & Dopitall993l) . 

2.3. Evolution of the Fermi Bubbles 

Some authors previously proposed, not a permanent CMZ 
wind, but a wind from a brief CMZ starburst 15 Myr ago; 
the shoc k front from t he wind could appear as radio/X -ray 
bubbles dSofuel l2000b lBland-Hawthorn& Cohenl HooH) . If 
the starburst wind just turned on, how long does it take to 
reach steady state and inflate the Bubbles? Suppose that the 
medium surrounding the CMZ is uniform, with density p ext 
and pressure P ext = p exl a 2 /2. The wind turns on, and starts 



sweeping up this medium at a shock. Once it reaches a ra- 



dius Rs e d « \/3M/(flv aC) p ex t), the mass of swept up material 
is greater than the mass in the wind, and the shock slows down 
- much like a supernova remnant does, except there is contin- 
uous mass injection. As long as \/3a < Voc, 7?sed < Ru f° r the 
conditions I have been assuming, /?sed ~ 1 kpc. 

After this point, I estimate the shock speed by assum- 
ing that the Bubbles are adiabatic, and that the kinetic en- 
ergy of the swept up material is a fraction k of the injected 
mechanical power. The radius of the shock is dR s /dt = 



6 K £//(pext^ 3 ), giving t = /? 5 /3[(3/50)p ext /(Ert)] 1 / 3 . It 
takes 



'inflate — 
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Myr (7) 



for the shock to reach the location of the termination shock 
(with pext = O.OOIoth cm" 3 and P e xt = 1000 K cm" 3 ), where 
the shock expands at a speed v s = 0.761k 1 ' 3 (v oo /ct) 1, ' 3 (T » 
230k 1 / 3 km s" 1 . The shock initially overshoots R, in this pic- 
ture, but then it oscillates back and forth before settling down 
at/?,. 

In my model, ?i n flate is longer than the minimum time H/E ps 
2 Myr for the starburst to do enough work to inflate the Bub- 
bles, where H = 7P ext [f2Z? 3 /3] is the enthalpy (c.f.. iCrockeri 
120121) . The reason is that the swept up material is not sim- 
ply pushed out of the Bubbles, but accelerated to v s , requir- 
ing more energy. On the other hand, my estimate o f H/E is 
~ 300 times lower than the 10 Gyr estimated by ICrockeri 
J2012I) . The major discrepancy is caused by three factors. 
First, the estimated volume of the Fermi Bu bbles in t h is dis- 
cussion is D,R^/3 «5x 10 66 cm 3 , whereas ICrockeri (12012b 
uses 2 x 10 67 cm 3 ; the latter va l ue is more accurate. Second 
and more importantly, Crocker] (12012b uses an external pres- 
sure of P/ks ~4x 10 4 K cm" 3 . This pressure is characteristic 
of the inner Milky Way disk, but the pressure in the halo (and 
most of the volume) could be much lower. Finally, ICrockeri 
(120121) uses the minimum outflow power of E ps 10 39 erg s" 1 , 
from the expected CR injection rate. In any case, it is clear 
that the CMZ can inflate the Fermi Bubbles. 

What if the wind suddenly shuts off? Then the Bubbles 
collapse at a speed ~ a, specifically the sound speed. The 
collaps e time is very long, tho ugh, alm ost 40 Myr. The large 
stellar dLaunhardt et al.l 120021) and gas (iMolinari et al.ll2011l) 
masses in the nuclear bulge are consistent with star-formation 
over several Gyr. Even if the output power of the CMZ fluc- 
tuates on timescal es of < 10 Myr, as extragalactic true nu- 
clear starbursts do dMavva et al.ll2004T) . the termination shock 
remains roughly in the same place. Thus, the Fermi Bub- 
bles may be a relative l y perm anent feature of our Galaxy, as 
iCrocker & Ahar onian (201 1) originally argued. 

3. NONTHERMAL EMISSION FROM THE SHOCK 

Because of the adiabatic losses, only a small fraction of CR 
energy density remains as the wind reaches its termination 
shock, regardless of any additional radiative losses. Yet nei- 
ther the CR nor the thermal energy disappears during adia- 
batic expansion; it simply converts into bulk kinetic energy. 
At the termination shock, that energy is converted back into 
random particle energy. Most of the energy goes into heating 
the gas, but an appreciable factor (> 10%) is in the form of 



relativistic particles accelerated at the shock. Thus, termina- 
tion shocks at the Bubble edges inject relativistic particles in 
situ, circumventing any losses. 

3.1. Estimate of '7- ray emission 

The mechanical power of the wind is just E. Suppose 
the shocks convert r\ = 30% of this power into CR protons, 
which are accelerated with a E~ 2 spectrum. I define \P = 
ln(£ , max /£ , m i n ), where E m [ n and E mix are the minimum and 
maximum energies of CR protons in the spectrum. Then the 
injection spectrum of CR protons is: 



> dQ P 

dE 



2.9xlO j5 ergs" 1 SFR ( , 1% . 3 *2(i, 



(8) 



where SFR .i = SFR/(0.1 M Q yr" 1 ), 770.3 = 77/O-3, and * 2 o = 
\&/20. The protons presumably diffuse out, or join a larger- 
scale Galactic wind, without radiating. 

Collisionless shocks also accelerate primary CR electrons, 
but > GeV electrons are thought to carry little of the energy. 
In the test particle approximation to diffusive shock acceler- 
ation, at energies > m p c 2 , the ratio of injected electrons to 
injected protons is 6 = {m p /m e ) {p - l) l 2 dBelllll978l) . Thus, for 
p = 2.0, we have 8 w 43: 

E 2 ^± w 6.7 x 10 36 erg s" 1 SFR .i 770.3*20- (9) 

dE 

Radiative losses for electrons are fast near the Galaxy, and 
the IC cooling time for an electron is 



tic = 3.1 Myr 



100 GeV 



U a 



1 eV cm 



(10) 



I assume that a fraction fie w 1 of the power goes into IC, 
with the rest mostly going into synchrotron. In the Thomson 
limit, the energy of an IC upscattered photon goes as E 1 oc E 2 , 
where E e is the relativistic electron energy. Thus, each dex 
of electron energy is stretched out into 2 dex of 7-ray energy, 
and we have E 2 dL/dE m (fic/2)E 2 dQ e /dE. The ave rage dis- 
tance of the Bubbles is 10 kpc, and they cover 0.8 sr (ISu et al.1 
120101) . So the flux from the Fermi Bubbles at Earth is 



, dNic 
dE 



220 eV cm" 2 sec" 1 sr" 1 /icSFRq. 1770.3*20- C 11 ) 



This m atches the 7-ray flux from the Fermi Bubbles (ISu et al.l 
l2OT0h . 

3.2. Magnetic fields 

Synchrotron microwaves from the Bubbles demonstrate 
that magnetic fields are present. Its 23 GHz bright- 
ness temperature is ~ 50 /jKelv at \b\ < 30° (iDoblerl 
l20T2at lAde et al. 120121) . From minimum energy arg uments 
dBeck & Krausd 20051) and 7-r ay to radio ratios (IDoblerl 
I2012bt iHooper & Slatveri 120131) . the magnetic fields are 
roughly 4-6 fiG. Th e radio "spurs" m ay have higher mag- 
netic fields, ~ 15 /jG dJones et al.ll2012l) . At higher la titudes, 
the microwave emissio n vanishes, indicating low B (Dobler 
127)123 lAde!tlini20!l . 

These magnetic fields pose two challenges to my model. 
First, the synchrotron losses cannot be much faster than IC, 
or else there is not enough wind power to explain the 7-ray 
emission. But for B < 6 /iG, fie > 0.5. 



Second, if the magnetic pressure is much greater than P ext , 
that could invalidate my calculations of R s . For B = 4 /iG, the 
magnetic pressure is 5000 K cm" 3 , with a similar pressure in 
CRs. But the wind ram pressure at 4 kpc from the CMZ is 
only 2000 K cm -3 . The tension is eased if the CMZ SFR was 
~ 0.2 M Q yr" 1 in the recent past. 

3.3. Ultra-high energy CRs from the shocks? 

The maximum CR energy is limited by the magnetic field 
in the shock, the size of the acceleration region, and the radia- 
tive losses the CRs experience. The sheer size of the Bubbles 
makes them excellent plac es for the accelerat ion of extremely 
high energy particles (c.f., Cheng et al. 2012). 

Electrons are highly radiative particles, so their maximum 
energy is set by IC or synchrotron losses. They gain a fraction 
~ Vqo/c of their energy per Larmor orbit, for an acce leration 
time of f acc w 20E/(3vooeB) that rises with energy (iGaisserl 
119901) . By setting f acc to the energy loss time from synchrotron 
and IC losses, I derive 



E r < 
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(12) 

U includes the sum of the CMB energy densities and the mag- 
netic field energy densities; Klein-Nishina effects suppress 
losses off starlight for PeV electrons, but electrons still radiate 
IC emission by upscattering the CMB. These 7-rays cascade 
down to TeV ene rgies due to 77 absorption within the Galaxy 
(Moskalenko et al. 2006). TeV and PeV electrons are a possi- 
ble source of TeV 7-rays from the Fermi Bubbles. 

For protons, the most important constraint on the maximum 
energy is that thei r Larmor rad ius cannot be larger than the ac- 
celeration region (lHillasll9 84). Requiring that E p < R max ZeB, 
I find 

£ ,< 5ZEeV (i^)(JL). ,13, 

It takes f acc w 4 Myr to reach these energies if B = 5 /iG 

Nuclei at these energies are subject to photopion and pho- 
todissociation losses off the far infrared radiation from the 
Milky Way. Photodissociation and photopionic 7-rays cas- 
cade down to ^100 TeV as they propagate to Earth. Photopi- 
ons can also decay into neutrinos, which could be detected at 
~ 10 PeV. 

4. DISCUSSION 

The starburst power source of this model is known, the 
wind's existence is fairly well established, and its termination 



shock is almost inevitable if the wind becomes supersonic. 
Because the injection sites are the termination shocks, there 
are no problems with cooling. IC emission does not require 
dense reservoirs of gas or long timescales, and is consistent 
with the wind's mechanical power. Finally, the shock could 
be a source of TeV IC emission, PeV neutrinos, and EeV CRs. 

There are several problems that must be addresse d. First, 
why does the e~ spec trum harden at low energies dSu et al.l 
120 lOt lAde et al.ll2012l) ? This could be a cooling break, but 
only if the electrons are just a few Myr old. Yet the CMZ 
is likely older (^> 10 Myr), so older electrons are probably 
present on some level. Have the older electrons simply dif- 
fused away? Perhaps the CMZ fluctuates on timescales of 
~ 10 Myr, but is steady on longer timescales. Then, the Fermi 
Bubbles' IC emission represents the electrons from the most 
recent pulse of star formation. Either way, low energy IC and 
radio emission should be visible surrounding the Fermi Bub- 
bles. Or is the break in the injection spectrum itself? 

Second, why are the Bubbles elongated? This indicates the 
termination shock radius depends on angle, which could hap- 
pen if the external pressure depends strongly on midplane dis- 
tance instead of distance to the Galactic Center, or if power is 
concentrated along the axis of the outflow cone. 

Third, thi s model cannot explain th e possible 7-ray jets dis- 
covered by ISu & Finkbeinerl (120 12b in the inner regions of 
the Bubbles, which could be from Sgr A* itself. The "jets" 
could instead be re lated to stellar clusters in the CMZ (c.f., 
iCarretti et al.H2.Q13b ; perhaps they are smaller cluster winds 
within the CMZ wind. 

Finally, the flat 7-ray surface brightness of the Bubbles is 
inconsi stent with IC e mission from a thin shell around the 
shocks dSu et alj|2010l) . If there is turbulence within the wind 
bubble, second order acceleration (perhaps of the adiabati- 
cally cooled initial CRs) might b e able to solve this prob- 
lem, but the model proposed by Mertsch & Sarkar (201 1 ) pro- 
poses that the shock is growing instead of steady-state; any 
turbulence must be generated before the shock itself. Alter- 
natively, if the wind fluctuates on Myr timescales, internal 
shoc ks may develop and a ccelerate CRs within the Bubbles 
(c.f.. lDorfi &~Br eitschwerdt 2012b. 
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